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Many  infectious  diseases  including  malaria  are  preventable,  yet  it  still 
poses  a  threat  and  thus  remains  endemic  in  many  communities  due  to 
lack  of  appropriate,  sufficient  and  timely  control  policies.  Strategies  for 
controlling  the  spread  of  infectious  disease  include  a  reduction  in  both 
infected  populations  via  treatment  and  a  possible  reduction  in  the 
susceptible  population  through  vaccination  or  sensitization.  In  this  paper, 
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we  gleaned  on  some  existing  model  and  thus  carry  out  a  modification  by 
incorporating  a  vector  reduction  parameter  as  a  new  control  strategy.  We 
determine  the  basic  reproduction  number  of  the  modified  model  and  also 
investigated  the  existence  and  stability  of  the  disease-free  equilibrium 
(DFE)  points.  We  showed  that  the  disease  free  equilibrium  state  is  locally 

asymptotically  stable  if  R0  <  1  and  unstable  if  otherwise.  This  shows 

that  if  R0  <  1 ,  malaria  can  be  controlled  in  the  population.  Finally, 

numerical  simulations  of  the  model  carried  out  using  the  fourth-order 
Runge-Kutta  numerical  scheme  in  Matlab  shows  that  malaria  can  be 
eliminated  in  the  shortest  possible  time 

1.  Introduction 

Malaria  is  an  infectious  disease  caused  by  protozoan  of  genius  Plasmodium  parasite  and 
transmitted  between  humans  through  bites  of  female  Anopheles  mosquitoes.  It  remains  one  of  the 
most  prevalent  and  lethal  human  infections  throughout  the  world.  An  estimated  40%  of  the  world's 
population  lives  in  malaria  endemic  areas  with  most  recorded  cases  and  deaths  occurring  in  sub- 
Saharan  Africa  [1],  Children  under  the  age  of  five  and  pregnant  women  are  the  most  vulnerable  to 
the  severe  forms  of  malaria.  Each  year  2-3  million  children  die  from  Plasmodium  falciparum 
malaria  and  up  to  500  million  people  throughout  the  world  suffer  from  malaria  clinical  disease  [2], 
Four  species  of  the  parasite,  namely:  P.  falciparum,  P.  vivax,  P.  ovale,  and  P.  malariae  infects 
human.  The  parasite  requires  two  hosts  to  complete  its  life  cycle  -  the  vector  female  Anopheles 
mosquito  and  human.  The  bites/bloodmeals  of  infected  mosquitoes  are  the  mode  of  transmission 
of  the  parasite  between  the  human  hosts.  The  symptoms  in  an  infected  human  include  bouts  of 
fever,  headache,  vomiting  flu-like,  anemia  (destroying  red  blood  cell)  and  malaria  can  kill  by 
clogging  the  capillaries  that  carry  blood  to  the  brain  (cerebral  malaria)  or  other  vital  organs  [3],  On 
the  average  the  incubation  period  of  Plasmodium  falciparum  is  about  12  days  in  humans.  Infection 
can  be  expressed  in  three  ways.  Prevalence  of  infection,  or  parasite  rate,  describes  the  proportion 
of  the  population  harboring  malaria  parasites.  Parasitaemia  describes  the  density  of  parasites 
within  a  host,  and  is  thought  to  be  an  important  factor  determining  the  severity  of  disease  in 
humans  [4],  Intensity  of  infection  describes  the  number  of  separate  infections  received  by  a  host  as 
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different  strains  of  parasites  differ  in  their  antigenic  properties  [5].  Thus,  intensity  is  certainly 
important  in  determining  the  level  of  acquired  immunity.  It  might  also  contribute  to  severity  of 
disease  [6],  possibly  by  determining  the  probability  of  a  human  becoming  infected  with  a  virulent 
parasite  strain. 

In  the  recent  years,  global  eradication  and  control  efforts  have  led  to  a  surge  of  activities  leading 
to  many  studies  and  publications  [7].  Control  strategies  and  intervention  program  have  been 
adopted  worldwide.  Some  of  which  include  the  use  of  anti-malaria  vaccines,  insecticides-treated 
bed  nets  (ITNs),  control  of  breeding  environment,  and  biological  control  among  others.  These  are 
largely  used  in  malaria  endemic  countries  especially  those  in  Sub-Saharan  Africa  and  have 
somewhat  led  to  the  reduction  in  the  spread  of  the  disease.  The  proportion  of  the  population 
sleeping  under  an  ITN  increased  from  less  than  2%  in  2000  to  an  estimated  55%  in  2015  (range: 
50-58%).  Ensuring  access  to  ITNs  has  been  critical  to  increasing  the  proportion  of  the  population 
sleeping  under  an  ITN.  Nearly  500  million  ITNs  were  delivered  to  countries  in  sub-Saharan  Africa 
between  2013  and  2015,  and  the  proportion  of  the  population  with  access  to  an  ITN  increased  to 
an  estimated  67%  in  2015  (range:  61-71%)  [3], 

1.1.  Role  of  Mathematical  Model 

Mathematical  models  play  a  key  role  in  the  control  of  malaria.  It  is  of  great  importance  and 
tasking  to  review  all  types  of  models  in  one  article.  In  this  article  a  historical  path  has  been 
considered,  and  an  attempt  is  made  to  take  into  account  some  of  those  mathematical  models, 
which  are  primarily  focused  on  the  transmission  dynamics  of  the  infection  in  the  host  and  vector 
populations,  using  the  epidemiological  compartment  modeling  approach  [8,  9].  [10]  derived  a 
model  where  humans  move  through  multiple  Susceptible  Exposed-Infectious-Recovered  (SEIR) 
stages,  where  a  history  is  kept  of  previous  infections.  They  include  a  sub  model  for  the  mosquito 
population  with  subdivisions  for  juveniles  and  adults.  They  used  the  steady  state  value  for  the 
adult  mosquito  population,  from  this  sub  model,  as  the  input  into  their  model  for  malaria 
transmission.  In  this  paper,  we  modified  the  existing  model  of  [1 1]  by  including  a  vector  reduction 
parameter  in  order  to  determine  its  impact  as  a  control  measure  for  the  spread  of  malaria. 

1.2.  Definition  of  Terms 

Susceptible:  The  number  of  individuals  who  can  be  infected  but  have  not  yet  contracted  the 
malaria  but  may  contract  it  when  exposed  to  its  mode  of  transmission 
Infected:  The  number  of  individuals  from  susceptible  compartment  who  have  been  infected 

of  malaria. 

Recovered:  The  number  of  individuals  from  infected  and  treated  compartment  who  have 

recovered  naturally  or  clinically  and  back  to  normal  status  of  health. 

Vaccination:  The  introduction  of  a  vaccine  or  serum  into  a  living  organism  to  confer  immunity 
WHO:  World  Health  Organization 

Ra :  The  expected  number  of  secondary  cases  produced  by  a  single  (typical)  infection  in 

a  completely  susceptible  population 

2.  The  Model  Equation 

We  present  here  the  model  equations  in  [1 1]  in  term  of  actual  population  and  their  parameters 


2.1. 

dSH 

dt 


Model  Equations  in  Proportion 

=  Knh  -  ahS"Iv  +  yRh  -  fihSH  -aSH 


(1) 
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H 

dt 

dRH 

dt 

dSv 

dt 

dly 

dt 


_  dbSHIy 

~  r^H  31 H  Mh^H 

(2) 

■  =  rIH  —  yRH  —  /uhRH  +  aSH 

(3) 

Cl  C'S  \;  I  I  ] 

=  KNV  MVSy 

^  H 

(4) 

acS  t/  I  tj 

(5) 

2.2.  Assumptions 

The  following  assumptions  are  made  in  order  to  formulate  the  equations  of  the  model: 

a.  The  infected  human  host  who  recovered  naturally  having  acquired  some  level  of  immunity 
are  moved  to  the  recovered  class. 

b.  Mosquito  never  recovers  from  infection,  as  it  is  regulated  by  mortality  of  its  individual. 

c.  Mosquitoes  bite  human  host  randomly  (independent  of  their  infective  status). 

We  incorporated  a  vector  reduction  parameter  k  in  the  vector  population.  Thus,  our  emerging 
model  equation  is  illustrated  as  follows: 


dS„  „  ,T  abSHIv  „ 

,  =KNh  m  +/Rh  Mi,Sh  oSh 

at  Nh 

(6) 

-C 

1 

1 

in 

t; 

i 

2  * 

Q 

II 

(7) 

dRH  r  ^ 

—  rlH  yRH  /j,hRH+aSH 
dt 

(8) 

dSv  acS^/1  tt 

-AVNV  Mv$v 

dt  Nh 

(9) 

<N 

1 

> 

1 

Etl 

2  * 

Q 

II 

(10) 

The  total  population  sizes  NH  and  Nv  can  be  determined  by  S  H  +  IH  +  RH  = 
Iv  =  Nv  from  the  differential  equations 

Nh  and  Sv  + 

1 

in 

1 

II 

(11) 

and 

~KNv  Sv{Mv+,Cl)  K  (/*v  +  ^2 ) 

dt 

(12) 

which  are  derived  by  adding  Equation  (6)  -  (8)  for  the  human  population  and  (9)  -  (10)  for  the 
mosquito  vector  population. 
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3.  Analysis  of  the  Model 

abSH  Iv 

In  the  model,  the  term  - denotes  the  rate  at  which  the  human  hosts  SH  get  infected  by 

N  H 

d  CS  y  I  zj 

infected  mosquitoes  Iv  and  - refers  to  the  rate  at  which  the  susceptible  mosquitoes  Sv  are 

N  H 

infected  by  infected  human  hosts  IH .  Since  it  is  easier  to  analyze  our  model  in  terms  of 

S  TJ 

proportions  of  quantities  instead  of  actual  proportions,  we  make  the  transformation  sh  = - , 

Nr 

IRS  I 

ih  =  ,  r,  =  — —  ,  5  =  — —  and  i  =  — —  in  the  classes  S„ ,  /„  ,  RH ,  Sv  and  Iv  in  the 


population  respectively  and  m  =  — — .  This  is  done  by  differentiating  the  fractions  with  respect  to 

Nr 

time  t  and  simplifying  as  follows: 


dsh  1  dSH  dNH 
dt  Nh  _  dt  h  dt 

=  \- abmsh iv  +  yrh  -  //,, sh  - ash  -  sh  [{Xh  -  pi h  ) -  Si b] 

=  Ah  —  abmshiv  +  yrh  -ash  -\sh  +Sshih 

=  K  I1  -  sh  )  -  abmshiv  +  yrh  -  ash  +  Sshih  (13) 


dih  1  dIH  .  dNH 
dt  Nh  _  dt  h  dt 

=  abmshiv  -  rih  -  5hih_pihih  -  ih  [{lh  -//*)-  Sih  ] 
=  abmshiv  —rih  -Shih_\ih  + 


=  abmshiv  -(r  +  S  +  Ah  )ih  +  d2h 

(14) 

Similarly, 

dr  /  x 

—  =  rih  ~\r  +  K  Fh  +  ash  +  oA 

dt 

(15) 

ds 

~r  =  (l  -  \K  ~  5v  (/*v  +k\  )]-  5v  V,cih  +  K  Uv  +  k2  )] 

dt 

(16) 

~r  =  Vlcih  +  K  (/',  +  K )]  +  K  [k  +  k2  X*'v  -  9  -  ] 

dt 

(17) 

subject  to  the  restriction  sh  +ih  +rh  =1  and  sv+iv  =  1 ,  we  note  that  the  population  sizes  N H  (t) 
and  Nv  (t)  do  not  appear  in  the  system.  Therefore,  using  the  relations  rh  =  1  -  sh  -  ih  and 
sv  =  1  -  iv  lead  to  studying  the  system  of  differential  equations. 

ds 

~  =  K^-sh)-abmshiv  +y(\-sh -ih)-ash  +Sshih  (18) 
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“77  =  abmshiv  -{r  +  5  +  Ah)ih  +  d2h 
at 

“T  =  (l  -  *v  X«C4  +  K  (/'v  +  *1  )]  +  *V  [U  +  k2  X'v  -  l)  -  ] 

at 


(19) 

(20) 


Now  we  intend  to  analyze  and  investigate  the  existence  and  stability  of  the  associated  equilibrium 
points.  Assuming  that  all  the  parameters  are  non-negative,  and  solving  for  the  equilibrium  points 
by  setting  the  right-hand  sides  of  Equation  (18)  -  (20)  to  zero,  the  system  takes  the  form  as  shown: 

A  l1  -  sh )- abmshiv  +  y{ 1  -  sh  ~  ih )- ash  +  5shih  =0  (21) 

abms*hi*v  ~(r  +  S  +  Ah  )i*h  +  Sf  =  0  (22) 

(l  -  iv  \acih  +  iv  (juv  +  /<, )]+  iv  [(//,  +  k2  \iv  - 1)  -  Av  ]  =  0  (23) 

In  the  absence  of  infection,  iv  =  ih  =  0 ,  so  that  Equation  (21)  yields 


\-\sh+r-Jsl-asl  =0 
K+r  =  s*h{K+Y+<x) 
s*  _  \±Y_ 

h  A h  +  y  +  a 


Hence  the  model  has  a  steady  state,  E0  called  the  disease-free  equilibrium  points,  where 

£,=  (/'l+r  .0. 0) 

Ah  +  y  +  a 

3.1.  Stability  of  the  Disease  Free  Equilibrium  Point 

To  establish  the  stability  of  this  equilibrium,  the  Jacobian  matrix  of  equations  (21)-(23)  is 
computed  and  evaluated  at  E0 . 

Let  Fx  =  Ah (l - s*h )- abmsJl  +  y(\  -  s*h  - ih )-as*h+  Ss*hi*h 

F2  =  abms*hi*v  ~{r  +  S  +  Ah  )i*h  +  d*2 

^3  =  (!  -  *v  lacih  +  K  (av  +  k\  )]  +  K  [(Av  +  k2  X'v  -  l)  -  Av  ] 

(  \  r)} 7  r -jp1 

— -  =  -Ah  -  abmi*v  -y-a  +  d*h  =  -\Ah  +  abmi*  +  y  +  a  -Si*h)  — -  =  -y  +  8s*h  — -  =  -abms*h 

dsh  dih  div 


OF ' 


=  abmi* 


dF, 


dih 

dF3 

dsh 

dih 

dly 

di. 


=  ac-aci„ 


OF: 


di 


A  =  ~(ach,  +Av)+{kl-  k2  )(l  -  2 iv* ) 


At  the  steady  states  of  the  model,  the  Jacobian  matrix  at  E  is  given  by 
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Jr  ~ 


{xh  +  abmi *  +  y  +  a  -  Si*h)  -y  +  Ss*h 
abmi*v  -QT  +  2  Si*h 
0  ac( 1  -  i* ) 


-  a  bins  l 
abmsl 


(aril  +  4  J+fo-  k2  )(l  -  2i* ) 


(24) 


where  QT  —  (r  +  S  +  Ah ) 


JF  = 


■(4+y  +  a)  -/  +  — - - 

4/,  +  7  +  a 


Evaluating  the  Jacobian  in  Equation  (24)  at  E0  gives 

A  +  7 
vA  +  Y  +  a 


■  abm 


0 

0 


Qt 

ac 


abm 

vA  +  Y  +  a , 
-4V  +{kx  -k2) 


(25) 


To  get  the  eigenvalues,  we  obtain  the  characteristic  equation 
Thus, 


■(4a  +  y  +  a)~  q  -y  + 


<^4 /,  +  b> 
4/,  +  7  +  a 


■  abm\ 


4  +r 

(  4ft  +  y  +  cc ) 


1 J  E0  9^1  — 

0 

Qt  Q 

abm 

f  4,+y  ^ 

1  O  I 

v4ft  +  Y  +  a. 

0 

ac 

-v 

+  (&j  -k2)~ ( 

=0 


(26) 


-  _(4ft  +  y  +  a)-< 


Qt 


ac 


abJ  K+T 

v4ft  +  y  +  a 

-4V  +  (^i  -k2)-q 


=0 


which  yields 

[-(4,,  +  y  +  a)-q] 


q2  +  [Cr  +  4V  -  -  &2  )]#  +  Qt  [4v  -  (^i  _  ^'2 )]  - 


2bcm{Xh  +  y) 
4a  +  y  +  cz 


=  0  (27) 


The  eigenvalues  of  the  characteristic  equation  are  given  by 


_(4ft  +  y  +  a), 
-  \_Qt  +  4V  -  (ki 


k2)\± 


[fir  +  4V  -(*i  "*2)]2  -4|er[4v  -(^  -*2)]-aW-^±0_ 
V  [  4,  +  y  +  a 


2 


i.e 


_(4ft  +y  +  a). 


~\Qt+  4V  -^2)]±dCr  +4y  -fc  ~^2)]2  ~4gr[4v  -fc -k2)ll-R0) 

2 
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where,  R0  = 


a2bmc 


r  -3  ,  a 

K  +  y 

+  y  +  a 


Qt  [Av  (kt  k2 )] 


,  is  called  the  basic  reproduction  number. 


The  eigenvalues  are  hereby  analysed 
i.e  qv  =  +  y  +  a)<  0 


~[Qt  +  K  -{ki  —  ^2 )]  +  Vl^r  +  *v  ~{ki  ~^)]2  ~4Qt[K  ~{k\  A) 

2 

-  [Qt  +Av~  (k\  -  k2 )]  -  -y/tCr  +  K  ~  ik l  “  ^2  F  ~  ^Qt  \Av  ~  (^i  “  ^2  W  “  ^  o  ) 


If  1  -  /?„  >  0 ,  then  Rn  <  1  and 


q  ;  -[Cr  +^v  -fe2)]  ,  V[Cr  +  4,  ~(*i  ~^2)]2  _0 


[Cr  +  _  (Xi  _  ^2 )]  Vfer  +  K  (^i  ^2 )]  F/o  ,  i  t?  i  Y1 

and^  <  i ^ - - - - - =  -[Qt+K  -(*i  -*2)J 

Therefore,  q2  <0  and  q3  <0  thus  establishing  ql  <0,  q2  < 0,  q3  <0 


Using  [12],  the  following  theorem  holds. 

3.1.  Theorem 

Giving  the  system  of  equation  in  (21)-(23)  and  that  r,  8,  Xh ,  y,  2V ,  kx ,  k2  >0  and  <1,  then  the 

disease-free  equilibrium  state  is  locally  and  asymptotically  stable. 

4.  Results  and  Discussion 


Table  1.  Parameter  values  and  references 


Exp  erim  ent 

V  allies 

References 

Parameters 

1 

2 

3 

4 

a 

0.05 

0.05 

0.05 

0.05 

[8] 

b 

0.5 

0.5 

0.5 

0.5 

[8] 

c 

0.05 

0.05 

0.05 

0.05 

[S] 

m 

0.10 

0.10 

0.10 

0.10 

[1] 

A 

0.02 

0.02 

0.02 

0.02 

[1] 

s 

0.33 

0.33 

0.33 

0.33 

m 

r 

0.3 

0.6 

0.9 

0.9 

Assumed 

A 

0.071 

0.071 

0.071 

0.071 

[8] 

Y 

0.140 

0.140 

0.140 

0.140 

[9] 

GC 

0.3 

0.6 

0.9 

0.9 

Assumed 

0.40 

0.40 

0.40 

0.40 

[1] 

0.24 

0.24 

0.24 

0.24 

[1] 

iv 

0.33 

0.33 

0.33 

0.33 

[S] 

&l’k2 

0.20. 

0.20. 

0.20. 

0.20, 

Assumed 

0.23 

0.23 

0.23 

0.23 
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Figure  1:  The  behaviour  of  the  model  when  r  =  0.3,  a  =  0.3  and  kx  =  0.20,  k2  =  0.23 

The  effect  of  implementing  treatment  and  vaccination  at  30%  is  shown  in  Figure  1  with  vector 
control  fixed  at  20%  and  23%.  While  treatment  and  vaccination  reduces  the  proportion  of 
susceptible  humans  we  see  that  the  disease  is  quickly  eradicated  from  the  population.  We  see 
clearly  in  Figure  2  that  increasing  treatment  and  vaccination  rates  to  60%  significantly  reduces  the 
susceptible  and  infected  humans.  And  a  further  increase  to  90%  improves  the  trend.  On  the  hand, 
vector  control  parameter  ki  and  k2  (kx  ^k2)  results  enhances  the  pace  at  which  the  vector 
population  is  eradicated.  In  particular,  when  K  <  k2 ,  the  infected  vector  will  continue  to  reduce  as 
the  time  progresses.  We  also  see  in  Figure  2  and  Figure  3  that  as  intervention  are  increasing 
deployed,  an  equilibrium  is  established  in  population  of  the  susceptible  humans  and  infected 
vectors  at  least  once. 


Figure  2:  The  behaviour  of  the  model  when  r  =  0.6,  a  =  0.6  and  kx  =  0.20 ,k2  =  0.23 
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But,  when  the  vector  control  parameter  is  suspended  as  indicated  in  Figure  4,  we  see  an  almost 
insignificant  impact  in  the  rate  at  which  the  infected  vector  is  reduced,  even  though,  the  treatment 
and  vaccination  rates  are  very  high. 


4.  Conclusion 

The  host  -  vector  malaria  disease  model  in  [11]  was  modified  by  including  a  vector  reduction 
parameter  for  the  dynamics  of  malaria  within  human  host  and  mosquito  vectors.  The  model  was 
analyzed  in  terms  of  proportion  of  quantities.  The  model’s  basic  reproduction  number  was 
determined.  The  stability  of  the  equilibrium  point  obtained  by  method  of  linearization  were 
analyzed  and  found  to  be  locally  asymptotically  stable  and  otherwise  if  1  -  R0  >  0 .  From  the 

result  of  our  simulations,  it  was  observed  that  treatment  and  vaccination  as  a  way  of  controlling 
the  menace  of  the  disease  in  the  human  population  is  quite  insufficient  but  with  the  combined 
effect  of  treatment,  vaccination  and  vector-population  reduction  as  an  additional  strategy,  malaria 
can  be  eliminated  in  the  shortest  possible  time. 


Nomenclature 

SH(t) 

5v(0 

JV(0 


a 

h 


r 

r 

S 


v 


A 

V 


Mv 


Number  of  susceptible  human  hosts  at  time  t 
Number  of  infected  human  hosts  at  time  t 
Number  of  partially  immune  human  host  at  time  t 
Number  of  susceptible  mosquito  vectors  at  time  t 
Number  of  infected  mosquito  vectors  at  time  t 
Number  of  female  mosquitoes  per  human  host 

Average  biting  rate  on  man  by  a  single  mosquito  (infection  rate) 

Proportion  of  bites  on  man  that  produce  an  infection 

Probability  that  a  mosquito  become  infectious 

Per  capita  rate  of  loss  of  immunity  in  human  hosts 

Rate  at  which  human  host  acquire  immunity 

Per  capita  death  rate  of  infected  human  hosts  due  to  the  disease 

Rate  of  recovery  of  human  host  from  the  disease 

Per  capita  natural  birth  rate  of  humans 

Per  capita  natural  birth  rate  of  the  mosquitoes 

Per  capita  natural  death  rate  of  the  humans 

Per  capita  natural  death  rate  of  the  mosquitoes 
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ki  Rate  at  which  susceptible  mosquitoes  are  killed 

k2  Rate  at  which  infected  mosquitoes  are  killed 


References 

[1]  WHO;  Global  Malaria  Programme:  Position  Statement  on  ITNs,  (2009). 

[2]  Engers,  H.D.  &  Godal,  T.  (1998 ).Malaria  Vaccine  Development ;  Current  Status. Parasitol, 
Today.Trend.Parasitol,  14,  pp:  56-64. 

[3]  World  Health  Organization  (2015).  "Q&A  on  artemisinin  resistance ".  WHO  malaria 
publications.  Archived  from  the  original  on  2014-07-20 

[4]  Ruiz,  D.,  Poveda,  G.,  Vlez,  I.D.(2006).  Modelling  entomological-climatic  interactions  of 
plasmodium  falciparum  malaria  transmission  in  two  Colombian  endemic -regions: 
contributions  to  a  national  malaria  early  warning  system.  Malaria  Journal ,  5  (66) 

[5]  Chiyaka,  C.,  Tchuenche,  J.  M.,  Garira,  W.,  Dube,  S.  (2008).  A  Mathematical  analysis  of 

the  effects  of  control  strategies  on  the  transmission  dynamics  of  Malaria.  Applied  Mathematics  and 
Computation ,  195,  pp:  641-662. 

[6]  Kawaguchi,  I.,  Sasaki,  A.,  Mogi,  M.  (2004).  Combining  zooprophylaxis  and  insecticide  spraying:  a  malaria- 
control  strategy  limiting  the  development  of  insecticide  resistance  in  vector  mosquitoes  Proc.  R.  Soc.  Lond. 
271,  pp:  301-309.  DOI  10.1098/rspb.2003.2575 

[7]  Alonso,  P.L  &  Teuscher,  T.  (2011).  Randomised  trial  of  efficacy  of  SPfi 66  vaccine  against  Plasmodium 
falciparum  malaria  in  children  in  southern  Tanzania ,  The  lancet,  pp:  1175-1181 

[8]  Aslan,  G.  &  Seyrek,  A.  (2007).  The  diagnosis  of  malaria  and  identification  of  plasmodium  species  hy 
polymerase  chain  reaction  in  turkey ,  pp:  87-102 

[9]  Barnwell,  J.W,  Collins,  W.E.  (2009).  'Plasmodium  knowlesi:  finally  being  recognized".  Journal  of  Infectious 
Diseases.  199  (8),  pp:  1107-8. 

[10]  Koella,  J.C.,  Boete,  C.  (2002).  A  theoretical  approach  to  predicting  the  success  of  genetic  manipulation  of 
malaria  mosquitoes  in  malaria  control  14(2),  pp:  56-64 

[11]  Adamu,  A.K.,  Wangercha,  WA.  (2013).  Modeling  the  spread  of  malaria,  IOSR  Journal  of 
Mathematics( IOSR-JM)e- ISSN : 2278-57 28 . Volume5 ,Issue5 (Mar. -  Apr .20 1 3), pp : 57-65  www.  iosr journals. org 

[12]  Van  den  Driessche,  P.  &  Watmough,  J.  (2002).  Reproduction  numbers  and  sub -threshold  endemic  equilibria 
for  compartmental  models  of  disease  transmission.  Mathematical  Biosciences ,  180,  pp:  29-48. 


112 


